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ABSTRACT 

In this paper, the validity of the Stokes-Einstein 
(SE) equation for estimating the transport 
properties (e.g. viscosity) of imidazolium-based 
ionic liquids [Camim]Br (n=2,4,5,6) was examined. 
To achieve this goal, the “hole theory” has been 
utilized for estimating the mean radius of the hole. 
According to the hole theory formalism, the mean 
radius of the hole is related to surface tension. To 
do so, the Fowler recipe of the Kirkwood-Buff 
(FKB) equation has implemented the calculation 
of Lennard-Jones's contribution to the surface 
tension of ILs. In the following, an assessment of 


the SE equation is depicted by drawing of = 


versus Z, Our test confirms the linearity of the diagram and then verifies the validity of the aforementioned approach. 


The results show that within the imidazolium family, as the cation alkyl chain length increases due to the increased 
Coulombic forces between the cation and the anion, surface tension decreases. 
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1. Introduction 

Ionic Liquids (ILs) as novel fluids have attracted great interest during the last few years [1]. Generally, ILs are a 
class of organic salt, which consists of organic cations and usually inorganic anions where one or both the ions are 
large, and the cation has a low degree of symmetry [2]. The physicochemical properties of the IL are mostly 
determined by the nature of the cation and or anion [3]. To optimize the use of ILs, especially in industrial applications, 
and design the new desirable ILs, a fundamental understanding of the physicochemical properties of ILs is very 
important [4]. ILs have extremely useful properties such as no detectable vapor pressure, wide liquidized range, high 
ionic conductivity [5], excellent thermal and chemical stability, and non-flammability [6]. 

ILs are being investigated for a variety of applications including extraction [7], separation [8], designable solvents 
for a variety of processes such as catalysis and organic synthesis, lithium batteries, fuel cells, super-capacitors, solar 
cells, material sciences, and engineering [9]. Thermophysical properties of fluids play fundamental roles in many 
diverse areas of science and engineering. They are divided into two categories: transport properties and equilibrium 
properties [10-15]. Indeed, the structural properties of ILs have received conspicuous attention [16-18] alongside both 
transport and equilibrium properties. 

Among the equilibrium properties ‘surface tension’ has great importance [11, 12]. Surface tension has the main 
role in distillation, absorption, separation, and extraction [19]. Knowledge of the surface tension properties and their 
relationship to the chemical structure of ILs [16-17] is essential to selecting the right IL, for a particular application 
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[20-21]. Despite the growing importance of ILs, the amounts of available correct experimental data on their transport 
properties of them are rare. The experimental data on surface tension for ILs reported in the literature are scarce and 
often conflicting; therefore, it is required to develop theoretical studies using correlation and predictive methods [7]. 
To do so, a microscopic understanding of the properties of ILs at a molecular level is necessary [22]. Therefore, 
we have a brief discussion on the theoretical models because they are simple and easy to use, and can greatly aid in 
predicting the properties and performance of ILs. Up to the authors' knowledge, there is no clear understanding 
regarding the prediction of the transport properties of ILs and different ionic structures by the SE equation. The hole 
theory (HT) is the most commonly used model for the investigation of ion movement in ILs at high temperatures [5]. 
Henceforth, the HT would need to consider the creation of holes within the fluid mixture instead. The aforesaid theory 
is one of the best models for the study of the motion of large charged molecules and low temperatures because the 
probability of movement of ions is smaller than in molten salts and also, an ionic association for charge transfer can 
be disregarded [23]. It is found that the simple lattice-HT theory can be usefully explained conductivity maxima [24]. 
The hole theory of liquids provides a simple and intuitive model which assumes that a liquid contains empty spaces 
or holes. In this model, space is considered to arise from thermally generated fluctuations in local density. The holes 
are randomly distributed and the radius of the average-sized hole (rH) is related to the surface tension (y). Once (rH) 
is known, the mobility of the hole can be obtained by the SE equation [25]. Therefore, the dire demand for 
thermophysical properties of ILs encouraged people to work out on these materials. Among the ILs, imidazolium- 
based ILs have received conspicuous attention. The selection of the imidazolium ring as a cation was due to its stability 
within oxidative and reductive conditions, the low viscosity of imidazolium ILs, and their ease of synthesis [26]. 
The surface tension of ILs is related to their pair potential function, radial distribution function, and liquid 
densities. We aim to evaluate a thermodynamic model for the prediction of the surface properties of imidazolium- 
based ILs by statistical mechanics. In this respect, contributions to surface tension from Lennard—Jones dispersion 
force were considered. Therefore, in this research, we evaluated the SE equation in combination with the hole theory 
for predicting the transport properties in particular surface tension of [amim]Br (ILs) (a= ethyl, butyl, pentyl, hexyl). 


2. Theory 

In this section, we have presented an argumentation about the main equation and model, which are used in this 
research. Previous research has shown that the SE equation is in good agreement with the diffusion behavior of ILs 
[5, 27]. It is becoming common practice to interpret the transport properties of a solute species using the SE equation 
and apply it to the investigation of the hole formalism in these materials [28-29]. The holes can be regarded in finite 
dilution and are independent of one another [25,30-33]. The SE equation is [31]: 

kT 
— 6nrn o 
Where D is the diffusion coefficient of the species, 7 is the viscosity coefficient, kg is the Boltzmann constant, T 


is temperature, and r is the effective hydrodynamic radius (EHR) of the species molecule [32]. According to Eq. 1, 
one can estimate the viscosity coefficient by knowing the diameter of the molecule. We postulated this diameter is 
given by the HT through EHR. 

It should also be mentioned that the primitive model was introduced [32] as an alternate view of the liquid phase 
and later extended by Bockris et al. [33], for use with molten salts. In this model, cavities occur at random positions 
and with random sizes because of thermal fluctuations, where the extent of cavity formation depends on the surface 
tension of the liquid [5]. In this theory, ions are assumed to migrate through empty spaces (holes) formed and 
destroyed by statistical fluctuations [5]. The HT is only suitable for ions movement that is smaller than the hole 
formed in its vicinity of it. The average size of holes in a liquid is calculated by the following equation [30]: 

kgT\2 
<r>=051 (=) (2) 

Where kg and T have their usual meaning and y is the surface tension of the liquid. 

The values of Lennard-Jones's (LJ) contribution to surface tensions of ILs are calculated through Fowler, 
Kirkwood & Buff (FKB) [12] via; 
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ij 3 alt o® r3 -d/o oè 
AA r | + gpa O (3) 
The reduced form of Eq.3 is as follows: 
TE o? a” F ped fa" Za í 
a 2r*2 4r*8 3r We 2r? (4) 


Values of surface tension via the FKB equation are utilized for estimating the mean radius of the hole according 
to the HT. Insertion of the mean radius of the hole to SE equation on the one hand, and the calculated values of 
diffusion coefficients (obtained from Mean Square Displacement (MSD)) on the other hand give the viscosity 


coefficients. All of the implemented quantities in Eq. 3, have been taken from Ref. [13]. Indeed, the dimensionless 
parameters and the LJ parameters used in this study can be defined as: 


kT 
T=— 5 
- (5) 
ee as 
n= (6) 
(me)2 
1 
Dm2 
D* = —; (7) 
O E2 
Ej = JEE (8) 
_atG; ý 
Oi; = a (9) 
Ae (10) 
yj 2 


The cross-interaction parameters between different ions were calculated using the Lorentz-Berthelot mixing rules. 
In Eq. 3, £ is the depth of the potential well, ø is the finite distance at which the inter-particle potential is zero, d is 
the diameter of the hard-sphere which relates to o of LJ fluid by the following equation [33]: 


di tri (74a) 


2 (11) 
1+¢c (ale) +1+¢3 (T/a) 


Oj 


3. Results and discussions 


The numerical values of coefficients of Eq.11 and ci—c3 are dependent on the type of RDF applied. In this study, 
c1=0.2977, c2=0.3316 and c3=0.001047 coefficients were considered. The LJ parameter of anion and cation of ILs 
were taken from Ref. [13]. The LJ's contribution to the surface tension of four imidazolium-based (ILs) [amim]Br 

2 
(a= ethyl, butyl, pentyl, hexyl) and the dimensionless surface tension (y* = -5 were calculated by the above- 


mentioned equations and parameters and are given in Table 1. 


Table 1. The Lennard Jones's contribution to surface tension of Cnmim]Br (n=2,4,5,6) ionic liquids and their reduced surface 
tension at 400 K. 


Tonic Liquids [Comim] Br [Cmim] Br [Csmim] Br [Comim] Br 
LJ 
Y * (mN/m) 6.953 3.642 4.840 2.341 
o*y 
r=) 0.206 0.203 0.396 0.173 
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As can be seen from Table 1, [C2mim]Br IL has the highest surface tension. The surface tension exhibits a 
decrease as the alkyl group is lengthened. In the case of reduced surface tension, the same trend is observed, except 
for [Csmim] Br. The larger the alkyl chain on the cation, the smaller the surface tension for the ILs as a result of the 
increased Coulombic forces between the cation and the anion. The other parameters such as temperature, type of 
anions, etc., also affect surface tension values. Also, It should be noted that to calculate the exact value of surface 
tension, one has to consider the other contributions from hard-sphere repulsion, and the electrostatic interactions to 
surface tension. The obtained values of < r > through Eq. 2, for ILs, are given in Table 2. 


Table 2. The obtained average hole radius of Camim]Br (n=2,4,5,6) ionic liquids. 


Ionic Liquid 1 ym 
r 
[Comim] Br 2.13 
[Camim] Br 1.43 
[Csmim] Br 1.77 
[Comim] Br 1.23 


As can be seen from Table 2, [Ci:mim]Br IL has the highest radius (the lowest value of ©). According to 


equation 2, hole radius and surface tension have an opposite relationship with each other. Therefore, as the chain 
length increases, the radius of the hole increases and the surface tension decreases. 
According to the best of our knowledge, there is no experimental data for the viscosity of ionic liquids of this 


. . Dx 1. ; we ais 
series. However, to evaluate the viscosity values, the plot of = versus 7 is depicted in Figure 1. 
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Figure 1. The plot of (Dxy)T) vs 1/r for Camim]Br (n=2,4,5,6) ionic liquids. 


As shown in Figure 1, the slope of (Dx7)T) vs 1/r for Camim]Br (n=2,4,5,6) ionic liquids is 7x10. According 
to the SE equation, if the slope of (Dx7)AT) vs 1/r is equal to kp/6z, the validity of the used values is confirmed. 


The plot of ae 


T* 


versus = is depicted in Figure 2. 
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Figure 2. The plot of (D*x7*)AT*) vs 1/r* for Cnmim]Br (n=2,4,5,6) ionic liquids. 
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Moreover, based on dimensionless equations, the slope of versus = is 1/6z and the slope of the reduced plot 


is 0.0531. The aforesaid results of the slope are adequate. 


4. Conclusion 
The tension at the surface of a liquid is one of the most notable intermolecular forces. Statistical thermodynamics 
provides a good opportunity to evaluate the thermodynamic model for representing the surface tension of ionic liquids. 


Tersely speaking, assessment of the SE equation is simply depicted by drawing of 2a versus - and aa 


* 


oO 
versus — , 
a 


simultaneously. Our test confirms the linearity of the diagram. Henceforth, the aforesaid approach for implementing 
the HT alongside the SE equation for ILs is fairly elaborate. Also, it is possible to calculate the viscosity of ILs using 
the Stokes-Einstein equation. The results indicated that the cation alkyl side chain length affects the surface tension 
of imidazolium-based ILs. In general, this work showed that knowing only the molecular parameters (e, ø), the surface 
properties of ILs can be evaluated. 


Appendix 
Expressions for surface tension 

Kirkwood and Buff's theory has explained the relationship between surface tension and intermolecular forces acting 
at an interface between two fluid phases based on statistical mechanics analysis. According to Fowler's approximation 
of KB theory, the surface tension of polyatomic fluids regardless of the contribution of the vapor density is given [12]: 


__p? fa aur) 
yar] G 


r 
In this Equation, r is the intermolecular distance, p is the density of liquids, u(r) is the intermolecular pair 


yg(r) r*dr (A. 1) 


potential, and g(r) is the radial distribution function (RDF). The common LJ pair interaction potential is: 


u(r) = 4ey a" 7 o (A.2) 


Where rj is the distance between particles i and j and both €, and oj are determined by using the Lorentz-Berthelot 
rules. The LJ potential parameters are used from Ref. [13] and the selected densities and diffusion of [Camim]Br 
(n=2,4,5,6) are taken from Ref. [34]. 


Authors’ contributions 
All authors contributed to data analysis, drafting and revising of the paper and agreed to be responsible for all the 
aspects of this work. 


Acknowledgments 
The authors thank the University of Neyshabur and the Shiraz University of Technology. 


Declaration of competing interest 
The authors declare no competing interest. 


Funding 


This paper received no external funding. 


Data availability 
Data will be made available on request. 


References 
[1] Z. Lei, B. Chen, Y.M. Koo, D.R. MacFarlane, Introduction: ionic liquids, Chem. Rev., 117 (2017) 6633-6635. 
[2] M.J. Earle, K.R. Seddon, Ionic liquids: Green solvents for the future, Pure Appl. Chem., 72 (2000) 1391-1398. 


Materials Chemistry Horizons | 2022, 1(2), 137-142 141 


Materials Chemistry Horizons 


[3] H. Ning, M.Q. Hou, Q.Q. Mei, Y.H. Liu, D.Z. Yang, B.X. Han, The physicochemical properties of some imidazolium-based 
ionic liquids and their binary mixtures, Sci. China Chem., 55 (2012) 1509-15018. 

[4] S. Zhang, N. Sun, X. He, X. Lu, X. Zhang, Physical properties of ionic liquids: database and evaluation, J. Phys. Chem. Ref. 
Data., 35 (2006) 1475-1517. 

[5] H. Zhao, Z.C. Liang, F. Li, An improved model for the conductivity of room-temperature ionic liquids based on hole theory, J. 
Mol. Liq.,149 (2009) 55-59. 

[6] D.K. Singh, B. Rathke, J. Kiefer, A. Materny, Molecular structure and interactions in the ionic liquid 1-ethyl-3- 
methylimidazolium trifluoromethanesulfonate, J. Phys. Chem. A., 120 (2016) 6274-6286. 

[7] H. Zhao, S. Xia, P. Ma, Use of ionic liquids as ‘green’ solvents for extractions, J. Chem. Technol. Biotechnol., 80 (2005) 1089- 
1096. 

[8] K. Friess, P. Izák, M. Kardszova, M. Pasichnyk, M. Lané, D. Nikolaeva, P. Luis, J. C. Jansen, A review on ionic liquid gas 
separation membranes, Membranes, 11 (2021) 97. 

[9] S. Zhang, K. Dokko, M. Watanabe, Porous ionic liquids: synthesis and application, Chemical Science 6 (2015) 3684-3691. 

[10] A. Eskandari Nasrabad, N. Mansoori Oghaz, B. Haghighi, Transport properties of Mie (14, 7) fluids: Molecular dynamics 
simulation and theory, J. Chem. Phys., 129 (2008) 024507/1- 024507/6. 

[11] M.M. Alavianmehr, S.M. Hosseini, B. Haghighi, J. Moghadasi, Surface thermodynamic properties of ionic liquids from new 
molecular thermodynamic model and ion-contribution equation of state, Chem. Eng. Sci., 122 (2015) 622-629. 

[12] A. Mulero, C. Galan, F. Cuadros, Fowler’s approximation for the surface tension and surface energy of Lennard-Jones fluids 
revisited, J. Phys. Condens. Matter., 15 (2003) 2285-2300. 

[13] M.M. Alavianmehr, S.M. Hosseini, J. Moghadasi, Densities of ionic liquids from ion contribution-based equation of state: 
electrolyte perturbation approach, J. Mol. Liq., 197 (2014) 287-294. 

[14] S.M. Hosseini, S. Aparicio, M.M. Alavianmehr, R. Khalifeh, On the volumetric properties of 2-hydroxy ethylammonium 
formate ionic liquid under high-pressures: measurement and molecular dynamics, J. Mol. Liq., 266 (2018) 751-761. 

[15] M.M. Alavianmehr, R. Ahmadi, N. Aguilar, M. El-Shaikh, S.M. Hosseini, S. Aparicio, Thermophysical and molecular 
modelling insights into glycerol + alcohol liquid mixtures, J. Mol. Liq., 297 (2020) 111811. 

[16] J.E. Enderby, G.W. Neilson, Structural properties of liquids, Adv. Phys., 29 (1980) 323-365. 

[17] O. Borodin, G.D. Smith, Structure and dynamics of N-methyl-N-propylpyrrolidnium bis(trifluoromethanesulfonyl)imide ionic 
liquid from moleular dyanmic simulations, J. Phys. Chem. B., 110 (2006) 11481-11490. 

[18] Th. Koddermann, D. Paschek, R. Ludwig, Molecular dynamic simulation of ionic liquids: a reliable description of structure, 
thermodynamics and dynamics, Chem. Phys.Chem., 8 (2007) 2464-2470. 

[19] P.J Carvalho, M.G. Freire, I. M. Marrucho, A.J. Queimada, J.A.P Coutinho, Surface tensions for the 1-alkyl-3- 
methylimidazolium bis (trifluoromethylsulfonyl) imide ionic liquids, J. Chem. Eng. Data 53 (2008) 1346-1350. 

[20] N.A. Ghani, N.A. Sairi, M.K. Aroua, Y. Alias, R. Yusoff, Density, surface tension, and viscosity of ionic liquids (1-ethy]-3- 
methylimidazolium diethylphosphate and 1, 3-dimethylimidazolium dimethylphosphate) aqueous ternary mixtures with MDEA, 
J. Chem. Eng. Data 59 (2014) 1737-1746. 

[21] X. Wang, Y. Chi, T. Mu, A review on the transport properties of ionic liquids, J. Mol. Liq., 193 (2014) 262-266. 

[22] H. Weingärtner, Understanding ionic liquids at the molecular level: facts, problems, and controversies, Angew. Chem., 47 
(2008) 654-670. 

[23] A.P. Abbott, R.C. Harris, K.S. Ryder, Application of hole theory to define ionic liquids by their transport properties, J. Phys. 
Chem. B., 111 (2007) 4910-4913. 

[24] C.E. Woodward, K.R. Harris, A lattice-hole theory for conductivity in ionic liquid mixtures: application to ionic liquid + water 
mixtures, Phys. Chem. Chem. Phys., 12 (2010) 1172-1176. 

[25] A.P Abbott, Application of hole theory to the viscosity of ionic and molecular liquids, ChemPhysChem 5 (2004) 1242-1246. 

[26] K. Ghandi, A review of ionic liquids, their limits and applications, Green Sustain. Chem., 4 (2014) 44-53. 

[27] P. Hapiot, C. Lagrost. Electrochemical reactivity in room-temperature ionic liquids, Chem. Rev., 108 (2008) 2238-2264. 

[28] M.A. Vorotyntsev, V.A. Zinovyeva, M. Picquet, Diffusional transport in ionic liquids: Stokes—Einstein relation or “sliding 
sphere” model? Ferrocene in imidazolium liquids, Electrochim. Acta 55 (2010) 5063-5070. 

[29] A.P, Abbott, Model for the conductivity of ionic liquids based on an infinite dilution of holes, Chem.Phys.Chem., 6 (2005) 
2502-2505. 

[30] J.O'M. Bockris, A.K.N. Reddy, Modern Electrochemistry, Vol.1, Plenum Press: New York, (1970). 

[31] T. Koddermann, R. Ludwig, D. Paschek, On the validity of Stokes—Einstein and Stokes—Einstein—Debye relations in ionic 
liquids and ionic-liquid mixtures, Chem.Phys.Chem., 9 (2008) 1851-1858. 

[32] C. D'Agostino, Hole theory as a prediction tool for Brownian diffusive motion in binary mixtures of liquids. RSC Adv., 7 
(2017) 51864-51869. 

[33] N. Mansoori Oghaz, B. Haghighi, M.M. Alavianmehr, E. Ghiamati, prediction of solubility parameters based on the explicit 
expression of statistical thermodynamics, J. Solution. Chem., 42 (2013) 544-554. 

[34] S. Yeganegi, V. Sokhanvaran, A. Soltanabadi, Study of thermodynamic properties of imidazolium-based ionic liquids and 

investigation of the alkyl chain length effect by molecular dynamics simulation, Mol. Simul., 39 (2013) 1070-1078. 


Materials Chemistry Horizons | 2022, 1(2), 137-142 142 


